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Calibration of Vibration Pickups by the 
Reciprocity Method 

Samuel Levy and Raymond R. Bouche 

The reciprocity theory for the relationship between mechanical force and velocity and 
electric current and voltage is presented for a linear electrodynamic vibration-pickup cali- 
brator having a driving coil, a velocity-sensing coil, and a mounting table. The theory 
takes account of flexibility in the calibrator structure and in the coils, electric coupling be- 
tween different parts of the same coil and between parts of one coil and parts of the other, 
and flexibility in the magnet structure. The theory shows what measurements are required 
in using a linear electrodynamic calibrator for the absolute calibration of vibration pickups. 

A description is given of the mechanical arrangement and electric circuitry used at the 
National Bureau of Standards in calibrating calibrators by the reciprocity method. A 
typical velocity-sensing-coil calibration curve is presented showing the effect of pickup mass 
on the calibration factor. Typical examples of measurement of calibration factor and me- 
chanical impedance of pickups are also presented. 

Practical limitations of the reciprocity method due to nonlinearity resulting from lack 
of tightness of mechanical joints, resonance effects, amplitude effects, etc., are discussed. 

1. Introduction 

An early application of the reciprocity method to the measurement of mechanical quanti- 
ties was that by R. K. Cook [1] 1 in 1940. He applied the method to the absolute calibration 
of microphones. This was followed in 1948 with papers by H. M. Trent [2] and A. London [3], 
who showed how the reciprocity method could be used for the absolute calibration of vibration 
pickups. The reciprocity method was applied to the calibration of electrodynamic transducers 
in 1948 by S. P. Thompson [4] and to the calibration of piezoelectric accelerometers in 1952 by 
M. Harrison, A. O. Sykes, and P. G. Marcotte [5]. 

Commercially available calibrators for vibration pickups have a mounting table that is 
connected by a relatively rigid internal structure to a driving coil and to a velocity-sensing coil. 
This internal structure is supported from the frame by flexure springs or guide wires. John ( ! 
Camm, formerly of NBS, in 1953 showed how the reciprocity method could be applied to the 
calibration of the velocity-sensing coil of such a calibrator in connection with work for the 
Office of Naval Research. In 1955 the authors of the present report, in connection with work 
for the Diamond Ordnance Fuze Laboratories, extended Camm's theory to take account of 
pickup mass and flexibility in the internal structure for the frequently occurring case where the 
driving coil, sensing coil, and mounting table are mechanically joined at a point. 

The present report gives a broader basis for the application of the reciprocity method to 
the calibration of the velocity-sensing coil of calibrators, by considering both the internal 
structure and the magnet to be flexible and by taking account of electric-coupling effects. 
It also shows how the calibrator, once its velocity-sensing coil has been calibrated by the reci- 
procity method, can be used to determine the calibration factor and to measure the mechanical 
impedance of vibration pickups. 

2. Reciprocity Theory for Flexible Vibration-Pickup Calibrators 

In the appendix the general reciprocity theory for a calibrator is presented. This theory 
is applicable to a calibrator operating in any frequency or amplitude range in which it can be 
considered to be a linear system. Neither the internal structure nor the magnet structure 
need be considered rigid. The theory is applicable, for example, not only at low frequencies 
and at frequencies near certain resonances, but also at frequeneies above axial resonance where 

1 Figures in brackets Indicate the literature references at the end of this paper. 
397105—56 5 227 



the driving coil and mounting table move in opposite directions. The theory is also applicable 
where there is electric coupling between different parts of the same coil and between parts of 
one coil and parts of the other. 

The positive terminals of the driving and sensing coils are designated as those having 
positive voltages when the mounting table is moving inward. Conversely, positive current 
in the coils produces outward velocity at the mounting table. 

When the calibrator is energized with current in the driving coil at any particular fre- 
quency, the calibration factor, F, is defined as the ratio of the voltage in volts in the velocity- 
sensing coil to the velocity in inches per second at the surface of the mounting table. It is 
shown in the appendix that 

F=a+bY p (1) 

where Y p is the pickup mechanical impedance (lb-sec/in.), and all the symbols represent 
complex numbers in the manner commonly used for alternating-current electrical theory. 
The calibrator constants, a and b, are determined by the following two experiments and com- 
putational procedure. 

Experiment 1: The equipment is shown in figure 1 (a). Attach weights to the mounting 
table and measure for each weight the transfer admittance in amperes per volt, G, between the 
driving and sensing coils: 

£=S (2) 

■^ 

where 

Io= current in driving coil. 

Eq = voltage generated in open-circuited sensing coil. (It is permissible to leave a volt- 
meter permanently connected to the sensing coil. If this is done, the open-circuit 
condition is obtained when no other connection is made to the sensing coil.) 
Experiment 2: The equipment is shown in figure 1 (b). Couple a second vibration exciter 
to the calibrator being calibrated at the mounting table, and measure the ratio, R, of open- 
circuited voltages generated in the sensing and driving coils: 

R=E S JE°. (3) 

Computational Procedure: Determine the ordinate intercept, J, and the slope, Q, of the 
function W/(G— G ) when plotted against the weight, W, attached to the mounting table in 
experiment 1, where G is the value of 6 when 1^=0. Constants a and b in eq (1) are then 
given by 

a=0.01711 -y/jccJR, b=GMlQTjR/(j<*7), (4) 

where co is the frequency in radians per second, andj is the unit imaginary vector. 

It is shown in the appendix that the vibration-pickup calibrator calibrated by the reci- 
procity method can be used to determine the mechanical impedance, Y p , of a pickup by meas- 
uring the admittance, G p , eq (2), with the pickup attached to the mounting table, and the 
admittance, G , when nothing is attached to the table, by use of the relationship 

X v ~ 3 V386 Al - Q(G P - Go)) W 

To carry out the calibration of a pickup the following procedure is used: The driving coil 
of the calibrator is energized at the desired frequency and the value of Y v is determined. The 
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Figure 1. Measurement of constants of vibration-pickup calibrator. 

(a) Weights, W 7 , are attached to mounting table of calibrator, C, in experiment 1. (b) External vibration exciter, E, is required to provide 
sinusoidal motion to calibrator, C, in experiment 2. 

calibrator is then driven at the desired amplitude as indicated by the voltage, Eg, in the sensing 
coil, and then the velocity is given by 

v=E*/F, (6) 

where F is obtained by substituting Y p in eq (1). The output of the pickup corresponding to 
this excitation is then measured. Ordinarily the term bY p in eq (1) for F is negligibly small, 
except at frequencies high enough to cause relative displacements between the mounting table 
and the sensing coil. Where it is known that bY p is negligibly small, the determination of Y v 
may be omitted from the calibration procedure. 
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3. Measurement of Constants 



3.1. Accuracy of Weights 

The weights attached to the table in carrying out experiment 1 at the Bureau are shown 
at W in figure 1 (a). The weights increase in 0.1 -lb steps from to 1 lb. They vary less than 
0.1 percent from their rated values. Their attachment surface has a stud that engages the 
mounting table and a contacting ring, }i in. wide, which provides a connection of high rigidity. 
A film of oil is wiped on this ring before engagement to eliminate air in the contact surface. 

3.2. Accuracy of Frequency Measuremenis 

The frequency was measured with a calibrated electronic frequency meter. The indicated 
audio oscillator frequency did not differ from the measured frequency by more than 0.2 percent. 

3.3. Measurement of Transfer Admittance 

The circuit shown in figure 2 (a) is used to measure the transfer admittance. With the 
switch in the "up" position, the values of r and r x are adjusted until the voltage drops \E l2 \ 
across terminals 1 and 2, and |£i 3 | across terminals 1 and 3, are equal as measured on a high- 
impedance voltmeter. The magnitude of G is then given by 



\G\- 



_ r J r r 1 +r 2 i 
rr 2 



(7) 



The accuracy of \G\ determined in this manner depends only on the accuracy of the circuit 
elements and the repeatability of the voltmeter reading, but not on the accuracy of the volt- 
meter. Values of the resistances are chosen to load the amplifier suitably. Typical values for 
this calibrator are r 2 , equal to 10 ohms, and r+r 1} approximately 1,000 times greater than r 2 . 

Voltage \E 2S \ across terminals 2 and 3 is then read from the voltmeter. The switch is 
lowered to the "down" position. The value of r x is adjusted until \E U \ equals l-E^I, and voltages 
|2? 13 | and \E^\ are read on the voltmeter. The phase angle, </> G , of the transfer admittance, G, 
is the phase angle of /£ with respect to E%. Because the currents in r, r 1 , and r 2 are in phase with 
Ioj phase angle <f> can be determined from a construction of a polygon of voltages as 

^cos-^l-lH}). (8) 

A suitable value of r 2 is used so that r^lO -6 is greater than 100; then the approximate value of 
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Figure 2. Circuits used in the measurement of constants 
of vibration-pickup calibrator. 

(a) Measurement of transfer admittance, (b) Measurement of voltage 
ratio. 
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<t> G can be determined from 

(<t> ) approx=90°-cos- 1 (l-| ||j). (9) 

The use of both eq (8) and (9) determines the magnitude and quadrant for <f> G . It happens, 
though, except near certain resonances, that cj> G is near ±90°. Small errors in \E A[i \ or \E 13 \ 
result in large changes in the angle 

when 4> (/ is near —90°. Difficulty in measuring the small value of \E 43 \ results in error when 
<t> G is near 90°. Therefore, eq (8) is used to determine the magnitude of <j> G and (9) its quadranl . 

3.4. Measurement of Voltage Ratio 

The circuit shown in figure 2(b) is used to measure the ratio of the open-circuit voltages of 
the driving and sensing coils when the calibrator being calibrated is driven by an external vibration 
exciter. With the switch in the "up" position and with r x set at 10,000 ohms, resistor r is 
adjusted until the voltages across terminals 1 and 2 and terminals 1 and 3 are equal. The 
magnitude of the voltage ratio, R, is then given by 

l^l=Ky)oo+-/ < 10 > 

(In the case of a calibrator for which 10,000 ohms is not effectively an infinite impedance across 
the driving coil, /'i in figure 2(b) should be increased to an adequate value and eq (10) corre- 
spondingly modified.) 

To determine the phase angle of B, \E n \ and \E 23 \ are measured. The switch is then 
put in the "down" position, r\ is adjusted until \E U \ equals |2?i3|, and |2? 43 | is measured. Then 
from a construction of a polygon of voltages 



022= COS 

and 



'H!i!> < u > 

*a=oos- l (l-| j^-ccs- 1 ^!). (12) 



It happens that <f> R is near either 0° or 180°. Equation (11) is insensitive for these angles, 
either because \E 23 \ is very small for <j> B near 0° or because the cosine is insensitive to small 
changes in \E 23 \ near 180°. Therefore, eq (12) is used to determine the magnitude of <t> R and 
eq (11) only its quadrant. 

4. Typical Results 

4.1. Calibration of Sensing Coil 

Results obtained in the calibration of the sensing coil of a typical calibrator having a 
nominal 50-lb driving-force rating are now presented. All measurements were made after 
thermal equilibrium was approached. The positive terminals of the driving and sensing coils 
were determined as those having positive voltage when the mounting table had inward velocity. 
Positive velocity is then outward. 

First, the transfer admittance, G, experiment 1, was measured for a sequence of weights, 
W, on the mounting table at each frequency for which a calibration was desired. Typical 
results at frequencies of 900 and 5,000 cps are presented in figure 3. In figure 4 are the cor- 
responding plots of W/(G— G ) against IV. The data were fitted by a weighted least-squnn^ 
procedure for a straight line. The weighting as determined by a consideration of the expected 
error in 1/(67— G ) was 1, 2, 3, . . ., 10 for W=0A, 0.2, 0.3, . . ., 1.0 lb. Intercepts J and 
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Figure 3. Variation of transfer admittance, G, with weight, W, on the mounting table for a typical vibration 
pickup calibrator at frequencies of (a) 900 cps and (b) 5,000 cps. 
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Figure 4. Plot of W/(G—G ) as a function of W for a typical vibration pickup calibrator at frequencies of 

(a) 900 cps and (b) 5,000 cps. 

slopes Q, determined by this procedure, are at 900 cps J=0. 1089 /— 94.2° lb-ohm, and Q = 
0.00767 /11.8° ohm, and at 5,000 cps, J r =0.01094 /-77.4° lb-ohm, and #=0.0076 5/100.4° ohm. 
The voltage ratio, R, when driven by an external vibration exciter, experiment 2, was then 
determined as R=0. 1090 /— 2.6° at 900 cps and i?=0.249 5/— 2.0° at 5,000 cps. Next the 
constants a and b in eq (4) were computed. Substituting their values in eq (1) gave calibration 
factors of the sensing coil with a pickup of mechanical impedance Y p on the mounting table, 
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Figure 5. Variation of calibration factor of 
velocity-sensing coil with frequency and pick- 
up weight for a typical vibration pickup cal- 
ibrator. 
O=o lb, x=H lb, # = i lb. 
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at 900 cps, ^=0.1402 /-3.4° +0.000674 /12.6° Y p v-sec/in., and at 5,000 cps, F=0. 1584/5.5° + 
0.001360 /93.1° F p v-sec/in. 

In figure 5 are shown the magnitude and phase angle of the calibration factor of the 
velocity-sensing coil of this typical vibration-pickup calibrator as a function of frequency for 
pickups having a mechanical impedance Y p corresponding to weights W of 0, 0.5, and 1.0 
lb (F p =jcoII7386). It is evident that, for this particular vibration-pickup calibrator, the cali- 
bration factor of the velocity-sensing coil is nearly independent of frequency and pickup 
weight up to 900 cps. The curve for 1 lb is dotted in the vicinity of 2,700 cps because the plot 
similar to that shown in figure 4 was not linear beyond 0.5 lb. The source of this nonlinearity 
could not be localized. This nonlinearity has not been found in other calibrators of the same 
make and model number calibrated at the same frequency. 

4.2. Calibration of a Variable-Resistance-Type Accelerometer and a Piezoelectric-Type 

Accelerometer 

Results obtained in the calibration of a variable-resistance-type accelerometer and piezo- 
electric-type accelerometer are now presented. In these calibrations the circuit shown in 
figure 2 (a) was used with the pickup terminals designated 6 and 7, in place of velocity-sensing- 
coil terminals 1 and 3. The method given in section 3.3 was then used for the measurement 
of G p , the value of G with the pickup attached to the mounting table. 

For the piezoelectric pickup, the value of G v obtained at 900 cps was G v =I^jEo = 
28.51 /91.8° mhos, and at 5,000 cps, G p =37. 9 /43.6° mhos. With terminal 6 grounded and 
terminal 7 replacing terminal 3, the same procedure gave, at 900 cps, I%/E p = 16.35 /2.5° mhos, 
and at 5,000 cps, I%/E P = 3. 201 /— 44.6° mhos, where E p is the output voltage of the pickup. 
From eq (5), the mechanical impedance, Y p , of the pickup was computed by using, at 900 
cps, g o =26.81 /92.0° mhos, and at 5,000 cps, g =16.51 /18.1° mhos, and the values of J and Q 
obtained in section 4.1. This gave, at 900 cps, Y p =2. 71 2/83. 3° lb-sec/in. and at 5,000 cps, 
Y p = 18.2 3/70.5° lb-sec/in. Using these values of Y v in eq (1), the velocity sensing-coil cali- 
bration factors at 900 cps and 5,000 cps, respectively, are 7^=0.1399 /— 2.7° v-sec/in. and F= 
0.135 7/9.2° v-sec/in. The pickup acceleration calibration factor, F p , is 

= (/?/fff) Fg 

where Fg/(jui) is the acceleration calibration factor of the sensing coil. Using this equation, 
at 900 cps, 7^=0.0166 /— 3.4° v/g, and at 5,000 cps, 7^ = 0.0197 /6.5° v/g. Figure 6 shows the 
complete calibration of the piezoelectric-type accelerometer. 

The calibration of the variable-resistance-type accelerometer, which was performed in a 
similar manner, is shown in figure 7. To compute the results on the variable-resistance 
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Figure 6. Calibration of piezoelectric-type accelerometer at a nominal acceleration level of 2 g. 
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Figure 7. Calibration of variable-resistance-type accelerometer at a nominal acceleration level of 2 g. 



pickup at frequencies other than those at which the vibration-pickup calibrator was calibrated, 
it was necessary to interpolate suitably for a, b, R, and Q, and compute J from eq 4. This 
interpolation was only used at frequencies below 900 cps. In this frequency range, such 
interpolation results in no appreciable error for this calibrator when the resonant frequencies 
of the flexures are avoided. The frequencies used for the pickup calibration were not those 
at which flexure resonance was present. 

5. Effects of Calibrator Construction on Accuracy- 
Inaccuracy in the reciprocity method, exclusive of electrical measurements, arises primar- 
ily from deviations from true linearity in the calibrator performance. The more important 
forms of nonHnearity are discussed below. 
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5.1. Resonance Effects 

Resonance accentuates inaccuracy due to nonlinearity because relatively small changes 
in the calibrator constants, which are caused by very small changes in the structure of the 
calibrator and small changes in the exciting frequency, can cause large changes in the amplitude 
of vibration. Even when such changes do not occur, the large amplitudes associated with reso- 
nance may exceed the linear range. The resonance of the driving coil relative to the velocity- 
sensing coil reduces the accuracy of the reciprocity calibration over a small frequency range 
because \R\ is very small and therefore difficult to measure, and because the large amplitude 
of the driver coil may exceed its linear range. Operation of the calibrator near resonance 
frequencies should therefore be avoided. Some common resonance conditions in their usual 
order of appearance with increasing frequency follow: 

a. Resonance of mounting table and coils as a rigid body on the guiding flexures. 

b. Local resonance in the flexures. 

c. Transverse vibration of the shaft connecting mounting table and driving coil. 

d. Longitudinal resonance of shaft connecting mounting table and driving coil. 

e. Local resonance in table or coils. 

5.2. Transverse Motion 

Transverse motion from any source invalidates the reciprocity method, which requires the 
mounting table to have uniaxial motion. Transverse motion occurs at the resonance frequencies 
of the flexures and at the transverse resonance frequencies of the shaft on which the moving parts 
are mounted. At resonance, the amplitude will vary substantially for small changes in fre- 
quency and the resonance may occur at a frequency that is a harmonic of the driving frequency. 
Flexures can be detuned by attaching small weights to them. 

5.3. Tightness of Mechanical loints 

The tightness of the mechanical joints in a calibrator will have lit lie effect on its perform- 
ance at low frequencies, where the moving assembly of mounting table and coils acts as a "rigid 
body." At frequencies approaching those for internal resonance in the assembly, looseness 
in the joints can affect behavior by, for example, changing the structural stiffness or by intro- 
ducing coulomb friction. At very high frequencies, joints normally thought to be tigbt may be 
a source of erratic behavior, which is difficult to eliminate as tightness will ordinarily vary with 
temperature and time. Coating all joints with oil improves their rigidity. 

5.4. Amplitude Effects 

The nonuniformity of the magnetic fields surrounding the coils is a primary source of 
nonlinearity. This effect tends to be greatest at low frequencies, where the displacement ampli- 
tude is large for a given acceleration. At high frequencies, the displacement amplitude is 
relatively small, and this effect is negligible. The change in equilibrium position of the moving 
assembly of a calibrator, when its orientation is changed, may cause a change in calibration 
factor due to change in the effective field strength. 

The structure of a calibrator will ordinarily experience small deformations in its linear range 
at permissible current intensities. At resonance, the deformations are larger. Little is known 
about the linearity of structural damping for a structure as complicated as a calibrator; how- 
ever, so long as the damping is small, its possible lack of linearity would affect only t he resonant 
response. 

5.5. Temperature Effects 

Changes in temperature cause moderate changes in elastic constants and in the tightness 
of joints, and in this way can affect the calibrator performance near resonance. Changes in 
temperature may also affect the field (lux density, and thus the electrical characteristics of the 
driving coil, but should not ordinarily affect the calibration factor of the velocity-sensing coil. 
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5.6. Purity of Electric-Power Sources 

If the magnets are excited by direct current containing an alternating-current ripple, this 
ripple may appear in the velocity output of the calibrator owing to currents included in the 
driving coil by transformer action from the field. Such a disturbance can readily be detected 
by exciting the field and observing the velocity of the mounting table with the driving coil 
short-circuited. 

Any harmonic content in the power supply to the driving coil will excite the calibrator at 
the harmonic frequencies, as well as at the primary frequency. The techniques described in 
this paper are applicable only for excitation at a single frequency. 

6. Conclusions 

It is concluded that the reciprocity methods described can be used for the accurate cali- 
bration of vibration pickups and calibrators having linear response. The accuracy is limited 
primarily by the accuracy of electrical measurement, by deviations from true linearity in con- 
struction, and by impurities in power supplies to the field and driving coils. The frequency 
range is limited theoretically only to frequencies for which the mechanical impedance of the 
weights attached to the mounting table can be computed. Ordinarily this limits the range 
to that at which the weights act as though they were rigid bodies. Practically, the calibrator 
used as an example in this report is most suited for use at frequencies of 900 cps and below, 
where its calibration factor is not affected by the mechanical impedance of most vibration 
pickups that would be attached to its mounting table. The calibrator is suited for use above 
900 cps, however, with some loss in precision primarily due to the electromechanical properties 
of the calibrator being less constant with respect to field-coil temperature than they are at 
lower frequencies. The upper acceleration range of most calibrators, when resonance is 
avoided, is below 50 g. This would be the limit for calibrating by the reciprocity method. 

7. Appendix. General Reciprocity Theory 

7.1. Reciprocity for Electric Circuit 

If a complicated electric circuit is considered to be made up of a number of meshes, the 
almost trivial fact that the common impedance between meshes i and k, for example, is the 
same as that between k and i is the basis for the proof of the reciprocity theorem. This proof 
is elegantly presented by Guillemin on page 152 of reference [6]. Guillemin states on page 276 
that the reciprocity theorem can be proved for both transient and steady-state performance. 
This means, for example, that if we impress a voltage in mesh 5 of a given network and measure 
the current, say in mesh 2, and then place the voltage in mesh 2 instead of in mesh 5 and meas- 
ure the current in mesh 5, we will find the current exactly the same in the two cases, both 
in magnitude and phase. 

7.2. Reciprocity for Mechanical System 

If a complicated mechanical system can be replaced by an equivalent system having 
discrete mass-points joined by springs and dashpots, the fact that the spring and dashpot 
connecting points i and k is the same as that connecting points k and i makes possible a proof 
of the reciprocity theorem for mechanical systems completely analogous with that for electric 
circuits. (Reference [7] gives a proof of this theorem for conservative systems.) In other 
words, this reciprocity theorem states, for example, that if we impress a force on point 5 of a 
given system and measure the velocity, say at point 2, we will find the same velocity in mag- 
nitude and phase when we place the force at point 2 instead of at point 5 and measure the 
velocity at point 5. It should be mentioned that it is possible to have a mechanical system 
consisting of discrete elements that is not reciprocal, for example, if it contains gyroscopic 
elements. 
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7.3. Reciprocity for Combined System 



In developing the reciprocity theorem for the combined electric-mechanical system of a 
calibrator, we will use matrix notation. This will permit us to consider the coils as flexible 
bodies and take account of electric coupling between coils, as well as deformation in the field 
and movable structure. The notation and theorems given by Frazer, Duncan, and Collar [8] 
will be used. 

Currents and voltages in the velocity-sensing coil will be indicated by S and in the driving 
coil by D. We will consider the coils subdivided into a sufficient number of segments so each 
moves as a rigid body and number them consecutively in the two coils. The subscript o will 
indicate terminal values. Complex notation used in alternating-current electrical theory will 
be used throughout. All symbols represent vector quantities, unless otherwise noted. An 
English system of units (inches, pounds, and seconds) will be used. 



Current /f at the terminals of the sensing coil is given by 



where 



E s E D 
I S o=$rs+fs"o+l2 s ]{E} 

^ oo <* ' oo 



(Al) 



Eo= voltage at terminals of coil S. 



Z^= voltage at terminals of coil D. 

Zoo S = impedance at terminals of coil S. 

Zof— transfer impedance from terminals of coil S to terminals of coil D. 






1 "1 where sf, is the transfer impedance from segment n to the termi- 



nals of coil S. 



{E} = 



E, 



E n 



E n is the voltage generated in segment n. 



A similar expression can be written for If. Both of these expressions can be expressed as 

{I } = [z ]{E } + [z]{E}, (Ala) 



where 



{/. 



U?J UfJ 



[«.]= 



[z]= 



1 1 
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where 



The current in any segment is given by 

{I}=[e]'{E.} + [Z\{E}, 

Ix" 
1% 

I n is the current in segment n, 



(A2) 
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2 mn is the transfer impedance from segment w to segment w, 



JL _L _i 

[z\= transpose of [2]. 

Velocity v at the"mounting|table, taken positive outward, is given by 

p.=|H-Zi[!/P]{I}, 

-*■ 00 

where 

i^ =the force at the mounting table. 
F 00 =the mechanical impedance at the mounting table. 
Z 1 =2.249X10~ 7 , a conversion factor. 
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In eq (A3), the product LiB n l n I n gives the force, in pounds, at segment n when I n is the current, 
in amperes. The velocity of the various segments is given by 



{v} = [y]'F +LAY\[k]{I}, 



(A4) 



where 



{v}- 



, v n is the velocity of segment n. 



[y]'= transpose of [?/]. 



[Y]= 



1 


1 


1 


2/n 


2/12 


2/1* 


1 


1 


1 


2/21 


2/22 


2/2^ 



_1_ J_ 

.2/ml 2/w2 



J. 

IJrnn — 



, 2/ TOn is the transfer mechanical impedance from segment m to 
segment 7i. 



The corresponding velocity in the magnetic field at the various segments is given by 

{v M }=L 1 [Y M ][k]{I}+[y"YF , (A4a) 

where 



{v M }-~ 



'vf 



v% 



v™ is the velocity of the magnetic field at segment n. 



pp 



yn yi2 
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'Jmn- 



, Vmn is the transfer mechanical impedance from the location of 
segment m to that of segment n in the magnetic field. 
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, yf n is the transfer mechanical impedance from the mounting table to the 
magnetic structure opposite segment n. 



The voltage generated in each segment is 

{E}= -Lmw+Lmm, 



(A5) 



where Z 2 = 2.540X10~ 8 , a conversion factor. In eq (A5), —L 2 B n l n (v n —v^) is the electromotive 
force generated in segment n, in volts, when the relative velocity between segment n and the 
magnetic field, in inches per second, is (v n —v„). Substituting (A2) into (A4) and (A4a), 
and substituting the results into (A5) gives 



where 




Letting 




where 






~1 




1 




1 



[1]= 



~FMk][y-y M Y-L l L 2 [K\[z]'{E }-L 1 L i [K\[Z]{E}, 
[K\ = [k][Y][k}-[k][Y M ][k]. 
[F- 1 ] = [l]+Z l i 2 [Z][Z] f 

<f 





, the unit matrix, 



(A6) 
(A7) 
(A8) 



.000... 1. 
[F _1 ]=inverse of [V]. 

{E} = -F L 2 [V}[k}[y-y M Y-L 1 L 2 mK\[z]'{E }. 
Substituting (A9) into (Ala), 

{I } = [z }{E }-F L 2 [z][V}[k][y-y M Y-L 1 L 2 [z}[V][K]M'{Eo}. 
Substituting (A9) in (A2) and the result in (A3) gives 

Vo^-InL i FAy-y M ][k][Z][V][k]\y-y M ] r +L 1 [y-y M ][k}[V]'[zY{E } 

J- oo 



using 



[VY=[l]-L,L 2 [Z][V][K\. 



(A9) 
(A10) 

(All) 
(A12) 
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The proof of eq (A12) is obtained as follows: By using (A8), 

[1}-L 1 L 2 [Z][V][K]=IK- 1 ]([K\-L 1 L 2 IKI[Z][V][K]) 

= [Z- 1 ]([iiq + (l-[F- 1 ])[F][Z]) (A13) 

= [K- l ][V][K}. 

Also, premultiplying both sides of (A8) by [K~ l ], 

[K-^V-^lK-^+L^Z]. (A14) 

As [K] is symmetric, [K~ l ] is symmetric. As [Z] is also symmetric, [K -1 ]^ -1 ] is symmetric. 
Its inverse [V][/£] is therefore also symmetric. The transpose of a symmetric matrix is itself or 

[V\[Kl=[K\'[V\'=[K\[V]'. (A15) 

Substituting (A15) in (A13) gives (A12). 

The current in the velocity-sensing coil is always zero, 

i?=[10]{J.}=0. (A1G) 

Substituting (A10) into (A16) and solving for Ogives 

F^zf^nvmiy-ry-E^^i-L^zf^wnK]^]') 

The current in the driving coil is 

/?=[0 1]{J,}. (A18) 

Substituting (A 17) into (AlO) and the result into (A18) gives 

/ z ss 

{%[z s W)[k][v-y M }'-[z D )[V][k][y-y»Y 

\ Z oo 

-L 1 L 2 e[2 s ][F]t/C]k1'[ Z D ][F]W[y-)/*]'-L 1 L 2 et2 5 ][F][A'][ 2 D ]'[^j[F]tfc][ ? /-- !/ »]') 
l-L,L 2 z?,f[z s ]F][tf][z s ]' 

(zk-j^r-L^ % [z*][V}[K][z s ]' + 2L l L 2 % [z"][V][K][z D ]' - L 1 L 2 [z D ][V][K][z"Y 

\ Z oo \ Z oo ) z oo Z oo 

+ mJUSfW [V] [K] ( [z s ] ' [z") - [z°] ' [z s ] ) [V][K] [z D ] ' ) 
+ E °~ l-LiL 2 z s J[z s ][V][K][z s ]' (A19) 

Substituting (A17) into (All) gives 

(± - UU\u ~ V*Wm\ V][k][y -y»]'- UU ^ [z*][V][K}[z s Y 

\ -* 00 J- 00 

+ LlLtz s Jlz s )[V][K][z s Y[y-y M ][Ic]m[V][k]ly-y»Y 
+ LiL 1 z°J[z*][V][k][ V -y><]'[z«][V][k][ y -y»Y) 



I?=F.L 2 - 



v =Fc 



1-L 1 L 1 z™[z8][V][K][zS]' 

(j£lz s ]lV][k][y-yi«Y-[z D ][V][k][y-yn' 

\ z oo 

+ L 1 L 2 z!f[z s ][V][K][z^[z D ][V][k][y 

' LiE °- l-LiLtzswmiKw (A20) 
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The coefficient of L 2 F in equation (A19) equals the coefficient of —L x Eo in (A20). Equa- 
tions (A19) and (A20) are subjected first to the condition v o =0 and then to the condition 7?=0. 
Eliminating E% in (A19), when v =0, and F in (A20), when i?=0, we obtain 



F 



un 



E D 

L 2 v 



7f = 



(A21) 



This is the reciprocity relation that exists in the combined electromechanical system. 
Equations (A20), (A19), and (A17) can be written for brevity as 

v =CF +AL l E° (A22) 

I D = -ALtFo+B^E? (A23) 

ES^HF.+NLtE?, (A24) 

where the constants A, B, C, H, and N depend only on the calibrator construction. 

In eq (1) a relationship is given for the calibration factor, denned &&F=E%/v , as a function 
of the mechanical impedance, Y PJ of the object on the mounting table. We will determine this 
equation in terms of the calibrator constants A, B, C, H, and N of eq (A22) to (A24). If the 
force, F , is only the reaction to driving a mechanical impedance Y at a velocity v , 

v =-F /Y. (A25) 

Substituting (A25) in (A22) and solving for F , 

A T 

^ =_jB? cq=r/r (A26) 

Substituting (A26) into (A24) and (A25) and forming the ratio Ef/v 9 , 

We see that N/A in eq (A27) is a in eq (1) and (NC/A)—H in eq (A27) is b in eq (1) when Y is 
Y P . 

a=N/A, b = (NCJA)-H. (A28) 

In experiment 1, section 2, the transfer admittance G=I%/Eo is determined for a series of 
weights, W, attached to the mounting table. Substituting (A26) into (A23) and (A24) and 
forming the ratio G, 

r _I? _ B+(BC+A*L t )T ,. . 

Ef N+(NC~AH)Y' { ' 

In experiment 1, Yis the mechanical impedance juW/g, where W is the weight on the mounting 
table, to the frequency in radians per second, and g the acceleration of gravity, 386 in. /sec 2 . 
Making this substitution in (A29) and forming W/(G—G„), where G is the value of G when 
F=0, we find 

W _386 7V 2 N 2 C-AHN w , . 

~G-G jo> A 2 NL 2 +ABH + A 2 NL 2 +ABH { j 
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We see that the intercept J and slope Q found in determining W(G—G ) experimentally are 

T 386 N 2 



ju A 2 NL 2 +ABH 



(A31) 



N>C-AHN (A , 

In experiment 2, section 2, the voltage ratio, R=E S \E®, is determined when I„ is zero. 
Eliminating F between (A23) and (A24) and setting J?=0, wc see that 

R=Ef/E»= HBLl+ .f NLlL2 (A33) 

Equation (4) is obtained when eq (A31) through (A33) are substituted into (A28) and L x 
and L 2 have tlie values given previously. 

From eq (A29), G =B/N, when Y is zero, that is, G is the value of G when nothing is 
attached to the mounting table of the calibrator. If G v is the value of eq (A29) when Y= Y p , 
where Y v is the mechanical impedance of a pickup to be determined, we find that 

v - iv {(x p —(jq) (A 34) 



~(A 2 L 2 N-+AHB)-(N 2 C-AHN)(G P -G ) 



The substitution of eq (A31) and (A32) into (A34) gives eq (5) for determining the mechanical 
impedance of a pickup attached to the mounting table of the calibrator. 



We are indebted to L. R. Sweetman for guidance in setting up the laboratory equipment, 
and to Richard Harwell, Jr., for the careful machining of the masses and the various fixtures 
used. Ruth Woolley performed the least-squares calculations and Carol Waldron prepared 
the figures. 
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